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"v^" ; Abstract 

O ' Experimental data obtained at the Yakutsk array after the modernization in 1993 are analyzed. The 

^ . characteristics of EAS longitudinal and radial development found from the charged particle flux and EAS 

! Cerenkov light registered at the Yakutsk complex array are presented. The energy spectrum of EAS obtained 

O ' from Cerenkov light and an estimate of the PCR mass composition are presented. 
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I. INTRODUCTION 



The measurements of charged particles (electrons and muons with > 1 • sec 6 GeV) and 
Cerenkov EAS radiation are carried out at the Yakutsk EAS array during more than 35 years. 
After the modernization in 1993 [HI], the Yakutsk array significantly improves the measurement 
accuracy of main EAS characteristics and increases a temp of statistics set in the energy range of 
10^^ — 5-10^^ eV. It has come about through the increase of the number of measurement stations 
and the decrease of separation between them. 

A spectrum of energy, dissipated by primary cosmic rays (PGR) in extensive air showers (EAS), 
have been obtained from these data and estimation of PGR mass composition was made 12, [j]. An- 
other method for estimation of primary particles mass composition will be presented hereinafter. 

There is a perspective method for analysis of primary cosmic rays (PGR) chemical composition 
based on conjoined analysis of longitudinal (cascade curve) and lateral (structural functions of 
electron, muon and Cerenkov components) development of EAS. With such a complex approach to 
measurement of shower characteristics here we have an opportunity of full reconstruction of PGR 
mass composition using specific mathematical techniques, for instance, inverse problem solving 
method [4;,^^, simplex method liTJ and so on. 



n. LATERAL DISTRIBUTION OF DIFFERENT EAS COMPONENTS 
A. Charged particles 

Fig. [T] shows the average lateral distribution functions (LDF) of charged particles for Eq ~ 
10^^ — 10^^ eV constructed by the method used for the Yakutsk array [8J. From fig. [T] it follows 
that at Eq ~ 10^^ — 10^^ eV LDFs are well measured in the distance interval of 15 — 400 m 
from a shower core, whereas at Eq > 10^^ eV they are only measured in the distance interval of 
50 — 1300 m. The curves are the approximation p{R) by the function ([T]) from the work [i9|]: 
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where A^s is the total number of charged particles at observation level, R^,^, is a mean square 
radius of LDF of charged particles. It is seen from fig. [3l |4]that the function ([T]) describes well 
experimental data both at Eq ~ 10^^ — 10^^ eV and in the showers at Eq > 10^^ eV. 

Thus, the function ([T]) can be used to determine in the showers at the highest energies where 
there are no measurements of particle density at small distances. 

On fig. [2] a comparison is presented between the lateral distributions of charged particles and 
muons with E'th > 1 GeV and calculation result from work [10]. In the work [10], a hybrid 
scheme of EAS simulation is used together with QGSJETOl model. As it follows from fig. [2l 
calculations give more slope function as for charged particles, so for muons at distances more than 
600 m from the shower axis. If calculated value of p(600) (QGSJETOl model) is used for shower 
energy estimation together with model-independent method for energy estimation used in Yakutsk 
experiment [[ill] then the energy estimated with QGSJETOl would be underestimated by the value 
equal to the difference of densities p(600) (see fig.O- 
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Figure 1 : Average LDFs of charged particles measured at the Yakutsk EAS complex array. 
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Figure 2: Lateral distributions of charged particles and muons. Circles are data from the Yakutsk experi- 
ment. Curves — from QGSJET model calculations for the showers generated by the primary proton. 



B. Cerenkov light 

Cerenkov light measurements at the Yakutsk array last for more than 35 years. Since the 
year 1993 there are 50 operating Cerenkov detectors with receiving area of photocathode 176 and 
530 cm^. Observation results for last years are presented on fig.[3l On the same figure calculations 
are shown from the work lll2ll (Dedenko et al.) for primary proton and zenith angle ^ = 0°. In 
the work a 5-level scheme for air shower generation is used. One can see a good correspondence 
in experimental data at medium and large distances from the axis. At distances < 80 m measured 
flux of Cerenkov light is less than one following from calculations. This discrepancy could be 
explained as with distinct mass composition so with different zenith angle. Experimental LDFs of 
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Cerenkov light are given for 6 = 17 — 18°. 

On fig. |4]our results (circles) are given in comparison with the data obtained by Haverah Park 
group [13]. The solid line is calculation from work [14J. There is a good correspondence with the 
experimental data. 
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Figure 3: The lateral distribution of EAS Cerenkov light. Curves — QGSJET (Dedenko et al.). 




ni. LONGITUDINAL DEVELOPMENT OF EAS 

Longitudinal development of high- and ultra-high energy air showers was reconstructed from 
Cerenkov light registered at the Yakutsk EAS array. For this purpose a method proposed in 
work [6] was used. Results of the reconstruction are shown on fig. [5] in comparison with cal- 
culations obtained for different hadron interaction models. Calculations have been performed for 
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primary protons and iron nuclei. Experimental data presented on fig.[5]are well described by mod- 
els with fast development like QGSJET in case of primary proton. It can be concluded that primary 
mass composition is a mixture of nuclei with varying portion of heavy nuclei. The most significant 
change is observed in the energy region of 10^^ — 10^^ eV (experimental data tend towards heavy 
composition) and above 3 ■ 10^^ eV where mass composition is closer to proton. 
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Figure 5: Average cascade curve of the shower development at E^q = 10 eV. Curves — different hadronic 
model calculations. 



IV. EMPIRICAL ESTIMATE OF EAS ENERGY AT THE YAKUTSK ARRAY 

The primary energy of a shower at the Yakutsk EAS array is calculated with the expression 

The energy scattered by electrons in the atmosphere above the observation level is given by the 
expression 

Sei = A:(X,P,)-F, 

Here, F is total flux of Cerenkov light from the EAS and k{X, P\) — is the coupling coefficient 
that represents the transparency of the real atmosphere and character of the longitudinal shower 
development, where Pa is spectral atmosphere transparency (SAT), calculated during lidar mea- 
surements. 

The energy of electrons at the observation level is calculated as 

E,, = 2.2 ■ 10'^ ■ iV,(Xo) ■ Aeff, 

where A^s(-^o) is the total number of charged particles at sea level and Aeff is the absorption mean 
free path of shower particles obtained from the correlation of the parameters Ns{Xo) and Q{R = 
400) at different zenith angles. Other components are: = e^- N^; E^i = (0.12 ± 0.09) ■ E^; 
= (5.6 ± 2.2) ■ E^; E, = (0.64 ± 0.18) ■ E^. 
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Figure 6: Estimation of the EAS energy by density ps(600), /3^(600) and (5(400). 



Using experimental data presented on fig.[6l the following expressions have been obtained at the 
Yakutsk array for energy estimation by density of charged particles and muons with i?th > 1 GeV: 



Ig^o = 17.68 + 0.98 ■lgp,(600) 
lg£;o = 18.32 + 1.12 ■lgp^(600) 



(2) 
(3) 



V. ENERGY SPECTRUM IN THE ENERGY REGION OF ~ 10^^ - lO^o EV 

In addition to charged particle surface detection, there is another technique used at the Yakutsk 
array — the air Cerenkov light measurement, which can be used to draw out the cosmic ray 
spectrum in independent way [fisll . The spectrum covers wide energy region 10^^ — 10^° eV. On 
fig- El fig- [8] Cerenkov spectra are compared to calculations with different models of cosmic ray 
propagation in the Universe [1^ 17]. It follows from fig.|7]and fig.[8]that galactic model describes 
well our spectrum in the region of 10^^ — 10^^ eV [,16.1 and metagalactic model [fivl] — above 
10^^ eV. As it is seen from fig. [8l in the region of 10^^ — 10^^ eV a boundary of transition from 
galactic cosmic rays to metagalactic possibly exists. This hypothesis requires further research. 



VI. MASS COMPOSITION OF PRIMARY COSMIC RAYS 

To interpret Yakutsk experimental data a set of X^ax and ps(600) calculated values obtained 
with CORSIKA simulation code (v.6.0, QGSJET model) was used. Calculations were carried for 
five primary nuclei (p. He, C, Si, Fe) at three primary energy values 10^^, 10^^, 10^^ eV [18]. In 
the work, two-dimensional probability densities F(Xmax, p(600)) were used with preliminary pro- 
cedure of standardization of the experimental data over the whole (Xmax, Ps(600)) data set. In the 
numerical implementation of this method, variables r and p were used instead of (Xmax, Ps(600)): 
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Figure 7: Spectra measured in Yakutsk using the air Cerenkov light. Curves are resuhs of anomalous 
diffusion model calculations \\m . 
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Figure 8: Calculated spectrum of extragalactic proton (solid curve) and of galactic iron spectra (dotted 
curves) 11711 compared to all-particle spectrum from Yakutsk array. 
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where — is a mean square error. 

For each of considered energy value and kind of primary nuclei (including nuclei joint in groups 
p + He, C, Si + Fe) probability distributions densities /(r, p) were constructed. The intersection 
of /(r, p) surfaces gives lines mi and m2 which optimally divide nuclei into (p + He), C and 
(Si + Fe) groups respectively. The simulations showed that with dividing of data into three groups, 
the effectiveness of nuclei group (p + He) to fall into the zone 1 and of (Si + Fe) to fall into the 
zone 3 is up to 90%. In the zone 2 a strong mixing between showers from different primaries 
occurs and a portion of carbon is ~ 50 % from all particles in this zone. 

The Xmax value characterizes a maximum of cascade development in individual shower and 
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p(600) is the density of particles at the observational level. 

Fig. m presents the results of multicomponent analysis of (X^ax, p(600)) data set obtained at 
the Yakutsk EAS array. It is seen that there is a correlation between the observed maximum of 
shower development and the density of charged particles 

The analysis was carried out for three values of energy, 2.4 ■ 10^'', 9.8 ■ 10^'' and 4.8 ■ 10^^ eV. 
In this notion, the point cloud represents standardized values, whose location regions characterize 
zones directly connected with the mass number of a primary particle. Lines represent borders of 
such zones. In this case, lines mi and m2 optimally divide nuclei into groups (p + He), C and 
(Si + Fe). One can see from fig. HI that the points are distributed over the zones non-uniformly. 




Figure 9: Standardized X^ax and peoo experimental data at different energy values, mi is a borderline 
between nuclear groups (p + He) and C, m2 is a borderline between nuclear groups C and (Si + Fe). 



Distribution of statistics over the energy intervals draws attention. The portion of (p + He) 
nuclei increases from 50 % to 53 % and a portion of carbon nuclei — from 23 % to 31 %. At the 
same time, the portion of heavy nuclei decreases from 27% to 16 % with growth of energy from 
2.4 ■ 10^^ eV to 4.8 ■ 10^^ eV. The error to recognize the nuclei for the energies of 10^'' — 10^^ eV 
does not exceed 30%. Such a distribution of nuclei in PGR does not contradict the conclusions 
about increase of the portion of protons and helium nuclei in the limit energy region made in our 
earlier works [19, 20l l2lll where other methods were used. 

On fig. [To] a portion of light, medium and heavy nuclei obtained in our analysis is shown. 
It is seen from the figure that at £"0 > 4 ■ 10^*^ eV portion of heavy nuclei decreases. From 
fig. [TT]it follows that in the energy region 10^^ — 10^^ eV (In A) value has its maximum value and 
after 5 ■ 10^'' eV starts decreasing. Such a change in mass composition in the energy region of 
10^^ — lO^'' eV might be caused by modernized diffusion mechanism of cosmic rays propagation 
in the Galaxy lfl6ll and possible influence upon the spectrum by lighter extragalactic component 
arriving from beyond the Galaxy IItIi . 
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Figure 10: Fraction of the different nuclei in the energy range Eq = 10^^ — 10^^ eV. 
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Figure 11: Mean mass numbers in cosmic rays of high- and ultra-high energies. Curves 
from anomalous diffusion model of cosmic rays propagation (Lagutin et al. 2004). 
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